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Abstract 

The likelihood of detecting individual discrete sources of cosmic rays depends 
on the mean separation between sources. The analysis here derives the minimum 
separation that makes it likely that the closest source is detectable. For super- 
GZK energies, detection is signal limited and magnetic fields should not matter. 
For sub-GZK energies, detection is background limited, and intergalactic magnetic 
fields enter the analysis through one adjustable parameter. Both super-GZK and 
sub-GZK results are presented for four different types of sources: steady isotropic 
sources, steady jet sources, isotropic bursts, and jet bursts. 

Keywords: cosmic ray; anisotropy; propagation; discrete source 

1 Introduction 

The rapid growth of the Pierre Auger Cosmic Ray Observatory's data set has prompted 
speculation that one or more discrete sources should soon be detected at ultra-high cosmic 
ray energies. Just as nearby stars are easily detected against the night sky, any nearby 
cosmic ray source should manifest itself as a cluster of cosmic ray arrival directions in the 
sky. At sufficiently high energies, magnetic deflections of charged particles from nearby 
sources should be small, so tight clustering of those arrival directions can be expected. 
At the highest energies, all of the observed cosmic rays should be coming from nearby 
sources because pion photoproduction (the GZK effect) eliminates high energy protons 
from distant sources. 

It is presently not possible to predict with certainty whether or not new observatories 
(e.g. Auger P, Telescope Array [2]) will identify and study individual discrete cosmic ray 
sources. The nature of the primary particles is still unknown and so is the mechanism that 
endows them with extremely high energy. Even if it is assumed that the highest energy 
cosmic rays are protons produced by discrete extragalactic sources, those sources might 
be distributed densely like galaxies or widely separated like the most powerful quasars. 
The mean separation between sources is crucial in evaluating the probability that the 
nearest source is detectable. 

Assumptions about intergalactic magnetic field properties are also critical in evaluating 
this probability. Proton trajectories are bent by magnetic fields, so the flux of particles 
from a discrete source is spread over a solid angle of the sky. The amount of blurring due 
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to intervening extragalactic fields depends on the distance to the source, on the typical 
magnetic field strength in intergalactic space, and on the effective coherence length for 
magnetic field directions. Properties of intergalactic magnetic fields are highly uncertain 
[31 m [5]. The better- known magnetic fields within the Galaxy will blur any extragalactic 
source somewhat, but at the highest energies the blurring of an extragalactic proton source 
due to the Galaxy's fields is less than the detector's point spread function (cf. Appendix 
B). 

Prospects for charged particle astronomy may be greatest above the GZK energy 
threshold for pion photoproduction [6]. Protons that are now detected well above the 
GZK energy threshold must have been produced within the last few hundred million 
years. Their sources are necessarily close enough, and their magnetic rigidities are great 
enough, to expect that their arrival directions must point close to the celestial positions 
of their sources. Moreover, the GZK effect eliminates the isotropic background which, at 
energies below the GZK threshold, has been accumulating for roughly 5 billion years. 

It is important to note that the favorable prospects for super-GZK charged particle 
astronomy stem from two independent considerations. It is a coincidence that magnetic 
blurring by galactic magnetic fields and by plausible extragalactic fields becomes unim- 
portant for protons approximately at the threshold energy where the GZK effect erases 
the isotropic background. Magnetic fields have nothing to do with pion photoproduction. 
If the GZK energy threshold were lower, charged particle astronomy would be favorable 
above the threshold even though intervening magnetic fields might spread any one source 
over a large solid angle on the sky, simply because the background would be gone. More- 
over, if magnetic fields were negligible below the GZK threshold, point source excesses 
would be detectable despite the isotropic background. With the two effects taking hold at 
approximately the same energy, it is assured that super-GZK astronomy is signal-limited, 
whereas sub-GZK charged particle astronomy is background-limited. 

Neither the GZK erasure of isotropic background nor the defeat of magnetic blurring 
by high rigidity has a sharp onset with energy. Instead of designating any particular 
energy as the dividing line between sub-GZK and super-GZK particle astronomy, this 
paper will rely on the characterization of super-GZK as being signal-limited and sub- 
GZK as background limited. Signal-limited analysis will be called super-GZK, whereas 
background-limited analysis will be called sub-GZK. 

The basic model to be used throughout most of this paper assumes that ultra-high 
energy cosmic rays are produced by a set of sources with mean separation R and that all 
of them have the same luminosity Q (cosmic rays per unit time). The density of cosmic 
rays in the universe today results from that emission rate density over an accumulation 
time T (which is limited by electron pair production below the GZK threshold energy 
and by pion photoproduction above the threshold). The individual source luminosity Q 
can be calculated from R, T, and the observed cosmic ray intensity /. For any assumed 
distance to the nearest such source, the expected flux is then determined. Its detectability 
at that distance depends simply on the cumulative detector exposure and, in the case of 
sub-GZK astronomy, how much solid angle is spanned by the magnetically blurred flux. 

It is easy to understand that the nearest source is more likely to be detectable if the 
sources are widely separated. Today's observed cosmic ray density is Q x T/R^ (the 
density of sources being 1/-R'^), so each source must have luminosity Q cc R^ in this 
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model. A larger mean separation R requires a greater luminosity Q per source to account 
for the observed cosmic ray density. The flux from a single source at distance r will be 
Q/inr^ oc R^/r^. Since the nearest source is likely to be at distance r R, the expected 
flux from the nearest source grows linearly with R. In the super-GZK (signal-limited) 
regime, this suffices to show that sources are more detectable if the mean separation of 
sources R is large. (This result pertains only when R is less than the GZK attenuation 
length.) For the sub-GZK case, one must worry about magnetic fields enlarging the solid 
angle over which the signal is seen by the detector. If the magnetic smearing is stochastic, 
that solid angle grows linearly with r (and hence R), but the noise (fiuctuation of the 
background within that solid angle) grows only as its square root. The signal-to-noise 
(detect ability criterion) is therefore enhanced by large values of R also in the sub-GZK 
regime. 

The analysis proceeds to find the minimum separation R that is needed in order for 
there to be greater than a 50-50 chance that the nearest source will be detectable. The 
answer depends especially on the assumed cosmic ray accumulation time T for particles 
above an adopted energy cut and on the number of arrival directions in the detector's 
data set that are above that energy cut. In the background-limited case, the answer also 
depends on a parameter that characterizes small-angle magnetic scattering in intergalactic 
space. For any set of assumed values for these parameters, one obtains the minimum mean 
separation of sources for which the detection of discrete sources is likely. 

The basic model assumes that the different sources of ultra-high energy cosmic rays 
are identical. They persist in time and emit cosmic rays isotropically. The model can 
be modified to study sources that emit in collimated jets, or transient sources that emit 
their cosmic rays isotropically in a burst, or bursts with jets. Part of the analysis can be 
performed with an arbitrary luminosity function rather than assuming that all sources 
have the same luniosity. 

No cosmic ray observatory has so far achieved full-sky coverage. Partial sky coverage 
is an obvious handicap in detecting discrete sources, for the detector may not get exposed 
to the brightest source. In addition, the analysis becomes more complicated without full- 
sky exposure. Section 10 indicates the necessary modifications for realistic calculations 
with partial sky coverage. The simplifying assumption in the bulk of this paper is that 
the detector has the same exposure to sources anywhere in a fraction of the sky denoted 
Jd- For a detector at one mid-latitude site, the effective fraction fo is approximately 1/2 
(cf. section 10). 

2 Notation 

E Cosmic ray energy, measured in EeV (1 EeV = 10^^e\^). In this paper, the energy 
E usually denotes the minimum energy used in the search for a discrete source. Due 
to the steep cosmic ray energy spectrum, most of the included cosmic rays are not 
much above that minimum energy. 

Q The luminosity of a cosmic ray source (number of cosmic rays per unit time). 

n The spatial density of cosmic ray sources (number of sources per unit volume). 
The meaning of n is somewhat different in the cases of bursting sources. (See the 
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description of t] in this section.) 
R The mean separation of cosmic ray sources {n = 1/R^). R is measured in Mpc 

(3.26 X 10^ light years), 
r The distance to a source, 
ro The distance to the nearest source, 
c The speed of hght. 

T The effective time (measured in years) over which the cosmic rays have accumulated. 
This depends on energy. For energies below the GZK threshold that are considered 
here, cT is limited by pair production to roughly 1500 Mpc. Well above the 
GZK threshold, cT is roughly 30 Mpc. The energy dependence of T through the 
threshold energy region is governed largely by the thermal spectrum of the CMB 
target photons [7]. 

/ Or /(> E). This is the integral intensity (cosmic rays per unit area per unit solid 
angle per unit time). "Integral intensity" means integrated over all energies above 
E. 

t The total operating time of the observatory, measured in years. 

Aq The ground area of the observatory measured in km'^. For a source at zenith angle 

6, the effective collecting area is A = Aq cos 6. 
£ The detector's exposure to a discrete source (and its part of the sky). This exposure 

is measured in km?yr. For example, a detector array with acceptance out to 60 

degrees zenith angle has an exposure to a source that passes through its zenith given 

by 

£ = detector jar ea x (live-time) x v^/27r. 

fn The fraction of the sky well exposed to the detector. For an observatory with 
uniform full-sky exposure, fo would be 1. For a detector at a single mid-latitude site 
with uniform acceptance to 60° zenith angle, the acceptance varies with declination, 
but an effective fraction //) = 1/2 is appropriate. See section 10. 

N The total number of arrival directions in a cosmic ray data set from a single obser- 
vatory site. For the case of a full-sky observatory with uniform celestial exposure, 
is 1/2 the total number of events in the data set. 

uj The solid angle over which the signal events arrive from a particular source. These 
may be expected to be distributed as a 2-dimensional Gaussian given by width a. In 
that case, the solid angle is given by u; = iira"^. 

M Noise, i.e. the amount of fluctuation in the expected background count. For high 
statistics, the noise is the square root of the background: 

M = sjh^. 

S Signal. This might be the number of counts in a target region above an expected 
background count. If a source is expected to produce a Gaussian distribution of 
arrival directions, then the signal S might be the Gaussian-weighted sum (minus the 
expected Gaussian- weighted sum from the isotropic background). 

S Signal-to-noise detection threshold. For example, if a 5-sigma detection is required, 
S would be set equal to 5. It is the number of sigmas deemed necessary to qualify 
for a positive detection. In a prescribed single trial, S = 3 might be appropriate [S]. 
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K The number of in-target arrival directions needed for a background-free positive 
source detection. For example, a cluster of 5 super-GZK arrival directions might 
constitute persuasive evidence for a source in a full-sky survey, so K might be set 
to 5. For a previously identified most promising candidate source, finding even 2 or 
3 new arrival directions in a small target region might be sufficient evidence. The 
number adopted for K depends on the search circumstances. 

K The angular diffusion coefficient that governs how the solid angle of magnetic blurring 
increases with source distance: uj = nr. 

B Extragalactic magnetic field strength measured in nanogauss. 

Z Electric charge of the cosmic rays, measured in units of the proton charge. 

p Larmor radius. To a good approximation, p = E/ZB, where B is the transverse field 
strength in nG, energy E is measured in EeV, Z is electric charge, and p is measured 
in Mpc. This formula is appropriate for estimations in intergalactic space. Within 
the Galaxy, the same formula can be used if p is measured in kpc and B is measured 
in /iG. See Appendix B. 

L The magnetic field coherence length measured in Mpc. This is a typical distance 
over which an intergalactic magnetic field affecting ultra-high energy cosmic rays can 
be regarded as having a consistent direction. 

/i The expected number of cosmic ray sources within a volume of specified radius. The 
radius is usually the distance limit at which any source is expected to be detectable. 

Qj The emission solid angle of a single jet (coUimated emission of cosmic rays). 

Q The solid angle on a sphere at distance r from a source over which the particles from 
a jet have spread. 

rj The spacetime density of bursting sources (rate of bursts per unit volume). For an 
accumulation time T, n = r]T is the spatial density of bursts that contributed to the 
present-day cosmic rays, and R = is the mean separation between the relic 

fossils of those bursts. 

W The number of cosmic rays emitted in a burst. For identical sources of bursts, this 

number is the same for all of them, 
r The time interval over which cosmic ray protons from an instantaneous burst arrive at 

Earth. The spread in time is caused by different trajectory lengths due to intervening 

magnetic fields. The time interval r is expected to increase with distance to the source 

and decrease with cosmic ray energy. 
a The coefficient that governs how the time spread r of received cosmic rays from a 

burst source increases with distance to the source: t = ar (cf. Appendix C). 
^ Used in Appendix A, ,^ is a dimensionless measure of distance given by ,^ = r/cT. 

3 The basic model 

The simplest model is that the observed cosmic ray intensity /(> E) is the result of 
isotropic emission by sources with identical luminosity Q and spatial density n over some 
history of T years. The density of cosmic rays is then given in terms of the sources by 
nQT. That same cosmic ray density is given in terms of the observed intensity as Anl/c. 
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Equating these two expressions yields the source luminosity 



^ = -7Kr2^- (2) 



where the density n has been expressed in terms of the mean source separation R such 
that n = The flux from a source at distance r is (5/47rr^, and multiplying this by 

the exposure £ gives the signal 

Exposure £ is the time-integrated perpendicular collecting area for flux coming from the 
source direction, and £ is measured in units of km?yr. 

Note: The flux Q/Anr'^ should be multiplied by e~~^l'^'^ if the effective accumulation 
time T in equation [1] is governed by a propagation attenuation such as GZK pion pho- 
toproduction. Including this factor prevents a simple analytic solution. The factor is 
irrelevant for the sub-GZK analysis. Omitting the factor leads to a valid result also in 
the super-GZK case in typical circumstances, as explained in Appendix A. 

The product 1£ is closely related to the total number of events in the data set. That 
number is a convenient parameter. Consider first a single observatory site such as Auger 
South, and let N be the number of events above some energy cut. To be specific, suppose 
this data set includes those events with zenith angles less than 60°. For a source that 
passes through the zenith of the detector, it is observable 1/3 of the operating time and 
the average collecting area is 

<A>=- I A^cos{d)de = ^Ao 

TT Jq 271 

where Aq is the full detector area and 6 is zenith angle. The exposure to the source is 
then £ = ^ < A >= ^Aq t, where t is the cumulative detector operating time. The 
number of events in the full data set is the cosmic ray intensity / times the product 
of aperture and operating time t. The detector's aperture (accepting events out to 60° 
zenith angle) is ^^o; so N = ^Aq I t. Combining this with the expression for £ yields 

2 /3 

I£ = N = 0.117 N. 

37r2 

For a full-sky observatory (a second site in the other hemisphere), the exposure to this 
source would be little changed, but the number of events in the data set would be twice 
as great. Using the approximation that a two-site observatory has uniform celestial ex- 
posure [1] , the calculation for the one source here would apply to all sources. The rule to 
be adopted, therefore, is that I£ = 0.117 A^, where A^ is the number of events in the data 
set for a single-site observatory or half the number of events in the data set of a full-sky 
observatory. (See section 10 for further discussion of non-uniform celestial exposure.) 

Super-GZK analysis: 
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Suppose K signal showers are needed to make a positive source detection in a signal- 
limited regime. Substituting K for S in equation [2] gives 



where tk is the distance at which the expected signal S would be equal to K. A source is 
detectable within that radius around us if it is in the fraction Jd of the celestial sphere. 
The expected number n of detectable sources is that exposed volume times the source 
density 

Substituting the previous expression for rx, this becomes 

^~ 3 ^cTK^ • ^' 

Denote by tq the distance to the nearest source. What condition on R ensures that tq is 
likely to be less than r/^-? The probability that ro is less than is the complement of 
the Poisson probability that there are sources within that volume when the expected 
number is /x: 

^^(^0 <rK) = 1- Poiss{0; /i) = 1 - e"''. 

Setting this final expression to be greater than 1/2 (i.e. the condition that detection is 
more likely than not), gives the condition fi > ln{2). Then substituting the above formula 
for yU yields this lower limit for the mean separation of sources: 

«>(|^P^f. (5) 

Alternatively, using IS = 0.117A^ and f£, = 1/2, the minimum separation Rmin for which 
the nearest source is likely to be detected is given by 



Rmin 4.1 ^ cT . 



(6) 



This equation quantifies the conditions for detectability in the signal-limited regime. 
The required minimum separation Rmin shrinks for any energy cut (fixed T) as the number 
of arrival directions increases (e.g. as an observatory's exposure increases over time). 

The expression for in equation [3] is derived without including GZK flux attenua- 
tion. The resulting expressions for /i and R are valid only if is much smaller than cT. 
As discussed in Appendix A, this is usually satisfied. The above expression for Rmin is 
therefore normally justified. There is also an upper limit on i?, since the nearest source 
should not be much farther away than the attenuation length cT. Appendix A shows how 
to find this upper limit Rmax for which detection of the nearest source is likely. 

Sub- GZK analysis: 
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The signal from a single source is presumed to be spread over a solid angle uj by virtue 
of small magnetic bends while passing through intergalactic space with some unknown 
spectrum of magnetic field strengths and randomly changing directions. By invoking 
the central limit theorem, or by analogy with multiple Coulomb scattering, it is natural 
to postulate that this magnetic blurring increases linearly with distance: lo = nr. This 
relation defines the diffusion coefficient k. (The RMS angle in any plane that contains 
the central direction to the source increases with the square root of the distance traveled. 
The solid angle increases linearly with distance because it is proportional to the product 
of two such angles.) 

The expected background within the solid angle u is luS. The noise is the fluctuation 
in this background. For large statistics, that is the square root of the background, so 



X = VlnrS. 

Using the expression for S above, the signal-to-noise ratio is 




K 



Now suppose a detection requires that S /N > S for some number S, and denote by 
the distance at which this signal-to-noise ratio occurs. The equation above gives 

A source is detectable within the radius around us. The expected number of sources 
within that radius is the volume times the source density: 

Denoting by tq the distance to the nearest source, as above, detection of the nearest source 
is likely if the probability is greater than 1/2 that ro is less than r^. That probability is 
given by the complement of the Poisson probability that there are sources within that 
volume when the expected number is /i. As above, this means /i > ln(2), which here 
reduces to 

R > {^f" {cT^f ^- (9) 

This is the constraint on the mean source separation R such that a detection of significance 
S is likely. 

Using I£ = 0.117A^ and = 1/2, the minimum source separation for which it is 
likely that the nearest source is detectable becomes 



Rmin = 1.4 (CT S)2 -. 



(10) 



The smearing by magnetic deflection is presumed to be the result of many small- 
angle scatterings as a charged particle passes through irregular magnetic fields between 
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the source and the Earth. The process is mathematically similar to multiple Coulomb 
scattering, and a Gaussian distribution of arrival directions can be expected. In that case, 
the solid angle uj is not a simple target region. Instead, each arrival direction should be 
weighted in proportion to the expected Gaussian distribution. The simple S/N analysis 
here still pertains, provided the weighting factor is taken to be Attct"^ times the Gaussian 
probability distribution. See Appendix D for further details. 

4 Fiducial calculations 

The inequalities [5], [9] and equations E], [10] give expressions for the minimum value of the 
mean source separation R such that the nearest source is likely to be detectable. The 
minimum mean separations depend on variables that are not known with certainty. Here 
some fiducial parameter values are adopted for illustration. Readers who favor other 
parameter values can readily scale the answers for their values. 

For the super-GZK estimate, the fiducial estimates here will be based on cT = 100 Mpc 
and = 50. There is some energy for which the accumulation time is (100 Mpc)/c, and 
any observatory should eventually detect 50 air showers above that energy. For it will 
be assumed here that there is reason to suspect the existence of a source at a particular 
location and that a single-trial test has been prescribed. A cluster of A' = 3 arrival 
directions at that celestial position would be a strong positive result if there are only 50 
events in a data set covering much of the sky. 

With these adopted values for the parameters cT, and A^, equation [6] gives Rmin ~ 
25 Mpc. The nearest source is likely to be detectable if the sources are separated by more 
than 25 Mpc on average. 

A fiducial calculation for the sub-GZK case requires adopting a value for k, the coef- 
ficient that controls how magnetic blurring increases with source distance. Wild guesses 
are allowed here as the properties of extragalactic magnetic fields are poorly determined 
by observations or theory [3] . One simplistic model is that the magnetic field has a typical 
strength B that is randomly oriented, but a particle experiences the same orientation for 
coherence length L. Its total path of length r is made up of r/L deflections from these 
randomly oriented fields. Appendix B shows that, in this simplistic model, k is given by 

K = ^ LZ'^ByE'. (11) 

If L = 1 Mpc, B =lnG, Z = 1, and E = 10 EeV, then k = ^ = 0.014 sr/Mpc. 

As an example, suppose cT = 1500 Mpc (a typical survival distance for nucleons 
against pair production at energies above the spectrum's ankle), and let A^ = 10^ arrival 
directions. As in the super-GZK fiducial estimate, suppose a prescribed test has been 
applied to a suspected celestial location, so S = 3 is statistically significant (a 3-sigma 
result for a single trial). Using cT = 1500 Mpc, N = 10,000, S = 3, and k = 0.014, 
equation dn] gives Rmin = 40 Mpc. The nearest source is expected to be detectable at this 
sub-GZK energy provided the sources are separated by at least 40 Mpc on average. 

The estimate for the detection distance rs would not be meaningful if the magnetic 
blurring at that distance (and for the energy cut E) were to produce a solid angle greater 
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than 27r. In this fiducial calculation with k = 0.014 sr /Mpc, the magnetic blurring for a 
source at 40 Mpc is 0.56 sr, corresponding to a circle of 24 degrees radius on the sky. 

The fiducial parameters adopted here may be quite wrong. They are presented only 
to illustrate the use of the formulas and provide explicit answers that can be easily scaled 
for different values of the parameters. 

5 Comments about the dependences 

In the super-GZK regime, the signal can be maximized by lowering the energy cut as 
much as possible without violating the signal-limited condition. For the actual mean 
source separation R, therefore, detecting the nearest source is made more probable by 
reducing the energy cut as much as possible in the signal-limited regime. 

In the sub-GZK regime, there may be little dependence of Rmin on the energy cut. 
At least for the simplistic model treated in Appendix B, k depends on energy as IZ-E^, 
which is approximately the same as the energy dependence of the intensity /(> E) or 
number of events A^. The factor -f^ is almost independent of energy, provided the detector's 
acceptance (aperture) is not growing with energy. Increasing the minimum energy reduces 
the magnetic blurring, but the effect on S/N is offset by the lower statistics. (There 
may be some advantage in raising the minimum energy in the analysis if the detector's 
acceptance does increase with energy.) Near the GZK threshold, energy dependence enters 
also through the accumulation time T. Raising the minimum energy in the analysis would 
then reduce Rmin- 

The statistical criterion for detection encoded in E enters quadratically in the inequal- 
ity [9] and equation [TOl As noted above, this number can be made relatively small if a 
careful sky survey determines a single best candidate source, its solid angle extent, and 
an optimal energy cut. A 3-sigma result for a single test with new data would then be a 
compelling result. 

Intergalactic magnetic fields might someday be studied using discrete sources of cosmic 
rays at known distances. For now, gross uncertainty is encoded in the single coefficient 
K, and any fiducial calculation based on an adhoc value of k should be treated with 
appropriate suspicion. 

Large detector exposure £ is obviously crucial for detecting discrete sources. Inequali- 
ties |5] and [9] show that the minimum value for the mean distance between sources in order 
for them to become detectable will shrink inversely as the exposure increases. A vast 
increase in exposure would lead to source detections, or else it would radically shrink the 
viable parameter space that describes the source distribution and magnetic fields. See 
section 10 for additional discussion about the dependence on exposure. 

The expected number of detectable sources jj, grows in proportion to <£^^/^ for the super- 
GZK case of signal-limited detection. A 10-fold increase in exposure results in a 30-fold 
increase in the number of detectable sources, and this is true even if that number is less 
than 1! Success in charged particle astronomy is all about achieving huge exposure. 

The 8'^/'^ dependence oi ji on 8 is closely related to the "logN-logS" relation of as- 
tronomy, in which the source count varies like flux~^^'^ for a homogeneous distribution 
of detectable sources in Euclidean space, independent of their luminosity distribution. 
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Since flux = K/£ in the analysis here, the number of detectable sources increases with 
exposure like 8^^"^. 

For background-limited (sub-GZK) analysis, the number of sources grows with ex- 
posure also, but it increases less rapidly with £. The number of detectable sources is 
proportional to 8^^^ . 



6 Sources with jets 

Jets are a common phenomenon among objects known to produce energetic particles. 
Relativistic bulk plasma motion is advantageous in accelerating particles to high energies, 
so it is reasonable to conjecture that the highest energy cosmic rays could be emitted in 
collimated jets. The previous assumption of isotropic emission from cosmic ray sources 
might be inappropriate. 

Suppose a single jet emits ultra-high energy cosmic rays in a solid angle VLj. The 
luminosity of every source is calculated as before: 

The measured signal from any one source is zero unless the observer is within the solid 
angle of the beam, which grows with distance r from the source by magnetic deflections. 
Let VL be the solid angle on the sphere of radius r which gets flux from the jet. Then the 
signal S is 

S = -9-£ = ^I^£ 
f2r^ Q cTr"^ 



Super- GZK analysis: 

A simple assumption is that the magnetic rigidity of super-GZK particles is high 
enough that Vt k, i.e. magnetic deflection does not significantly increase the solid 
angle of the jets in transit between the source and Earth. The signal S at distance r is 
then given by the previous equation with f2 replaced by Vtj. 

For a given exposure £ and cosmic ray intensity / accumulated over T years, K 
particles would be expected at distance tk'- 

The probability is r2j/47r that Earth is in the beam of any one source, so the density of 
viewable sources is 4^^. The expected number of detectable sources (i.e. pointing at us) 
within the exposed volume of radius tk around us is this: 

^~ 3 i?3-3v^^ cTK ^ ■ ^^^^ 
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As in the basic model above, the probabihty that the nearest source is closer than 
tk is greater than 1/2 provided /i > ln2. Using the above expression for /i converts this 
inequality to a lower limit on the mean separation of sources: 



Detection of the nearest source is more likely than not provided the sources have at least 
this mean separation. Notice that the minimum separation is decreased relative to the 
basic model by the cube root of the jet opening solid angle (as a fraction of the isotropic 
An solid angle). 

A fiducial value of Qj will here be taken to be 0.01 sr, corresponding to a jet opening 
cone of 3.2° half-angle. Together with the previously adopted fiducial values (cT = 
100 Mpc, N = 50, K = 3), the fiducial calculation here gives Rmin = 2.3 Mpc, suggesting 
that the nearest source should be detectable even if the sources are distributed more 
densely than normal galaxies. 

The minimum separation Rmin is substantially smaller than in the basic model because 
the signal is strong when looking into a jet. The total luminosity of every source can be 
relatively weak (as must be the case if the density of sources is high), and an individual 
source can nevertheless be detected far away. (For Rmin = 2.3 Mpc in this fiducial 
estimate, the maximum detection distance is 17 Mpc, but that detection distance 
grows in proportion to R^/"^). 

In the super-GZK regime, jet sources are easier to detect than isotropic sources if the 
mean separation between sources is small, but they are more difficult to detect if the mean 
source separation is large. Since tk/ R is larger for jet sources than for isotropic sources, 
it may happen that there are many sources within the GZK volume but none of them 
is pointing at us. For example, suppose R = 25 Mpc (the minimum mean separation 
for detectability in the fiducial estimate for the isotropic case). The expected number 
of sources within 100 Mpc of Earth in the exposed half of the sky would be 134. For 
Qj = 0.01 sr, the probability is only = 8 x 10~^ for Earth to be in the beam of any 
one source. Although there are 134 sources readily detectable if any points at us, the 
expected number pointing at us is only 0.1. Detection is unlikely for a large mean source 
separation. 

Sub- GZK analysis: 

Although super-GZK particles can be assumed to maintain their directions in transit 
from source to Earth, that is not expected in the sub-GZK regime. As sub-GZK particles 
get farther from the source, the effective angular extent of the jet is dominated by their 
magnetic deflections rather than the emission angle Qj of the jet itself. Here it is therefore 
assumed that, at the distance to Earth, the flux from a single jet is spread over a solid 




(14) 



Using IS = 0.1 17A^ and //) = 1/2 as before, this becomes 




(15) 
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angle Q = Snr. (See Appendix B for an explanation of the factor 3.) The signal at 
distance r becomes (cf. equation [2]) 

Att is 



3 cTk, 

The background does not care that the luminosity of individual sources is in jets. It is 
given, as in the isotropic case, by luS. Its square root gives the estimated fluctuation in 
the background. Substituting nr for a;, that noise is 



A/" = ViKrS. 

The signal-to-noise is then 



3 V cTr^l'^ 

The analysis proceeds as in the case of isotropic sources. Adopting a value S for iS/A/" 
gives an expression for r^: 

This is the radius at which a source should be detectable with signal-to- noise equal to S. 

The expected number of detectable sources (beamed at us from the exposed fraction 
of sky within the radius is 

^r'^ndr = fo J ?>Kr'—dr = fj,-Kri/R\ 

Using the foregoing expression for r^, the expected number of sources closer than 
becomes 

4 3 K CI h 

As in the case of isotropic sources, the probability that the nearest source is detectable 
(closer than r-^) is given by /i > ln{2), which here reduces to 

R > (iM)7/3(A)8/3(^)8/3^5/3_ ^g) 



3/d ' Mvr' 

Substituting 0.117A^ for IS and 1/2 for fc as previously, this becomes 



Rrmn = 1.6(^)^/^(cT)«/3«:^/^ 



(19) 



Using again the fiducial parameter values S = 3, cT = 1500 Mpc, N = 10, 000, k = 
0.014 sr/Mpc for sub-GZK analysis, this formula gives Rmin = 33 Mpc. With this 
minimum source separation, sources are detectable (S > 3) out to = 47 Mpc, and the 
arrival directions are spread over solid angle uj = 0.66 sr (26° cone half-angle). 

The sub-GZK regime offers detection capability for jetted sources with large mean sep- 
aration even though super-GZK detection might be unlikely in that case. The two regimes 
are complementary. Super-GZK astronomy is likely for densely distributed sources with 
jets, whereas sub-GZK astronomy is likely for more widely distributed sources with jets. 
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7 Isotropic burst sources 



It is unlikely that the sources of high energy cosmic rays are permanent. The universe 
is dynamic, and even the most powerful active galactic nuclei may be temporary feeding 
episodes of supermassive black holes. For temporary sources, the source density n should 
be interpreted as the mean density of active sources at any time. The formulas of the 
preceding sections should then be appropriate. 

Some modifications are needed, however, if cosmic rays are produced in brief bursts. 
Due to the magnetic wandering of particles en route to the Earth, the duration of flux 
from an instantaneous burst increases with distance from the source. The duration r (at 
any cosmic ray energy) should increase linearly with distance: t = ar (where a is energy 
dependent, c/. Appendix C). In this paper, cosmic ray sources are regarded as "bursts" 
if their emission lifetimes are not long compared to the time spread r expected at Earth 
for a source that emits all of its cosmic rays instantaneously. Let W be the magnitude of 
a burst, i.e. the total number of emitted cosmic rays. If an isotropic burst is at distance r 
from the Earth, the flux over time r is W/ (47rr^r). Substituting ar for r and multiplying 
by the exposure £ gives the signal, 

s 



4:7iar^ 



Let r] denote the spatial density of bursts per unit time. If T is the accumulation time 
for cosmic rays, then n = rjT is the fossil density of bursts that have contributed to the 
accumulated cosmic ray density, so 

Air 

—I = nW = r]TW. 

c 

Solving this for W and inserting the result into the expression for the signal S gives 

s 



ac rjT r^ 

An estimate for a derived in Appendix C is 



a 



36cE2 ■ 
Super- GZK analysis: 

As usual, suppose K events are needed for a detection: 

S = K ^ rl= (20) 
ac rjl K 

The density of burst sources at distance r with flux now "on" is rjr = rjar. The number 
of detectable bursts closer than rx is therefore 

rrK 

/i = /x) / Airr^ rj ar dr = TrfnCtrfr'^. 
Jo 
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Using the foregoing expression for and = l/irjT), this can be written as 

= ^ )^/^ R- (21) 

ci ac K 

The probabihty P{rQ < tk), that the nearest (on) source is closer than r^, is greater 
than 1/2 (detection hkely) provided 

/i>M2) ^ i?>McT(ac)^/3(^)^/^ (22) 

Substituting 0.117A^ for I£ and 1/2 for as before, this becomes 

(23) 



Rrmn='^-^{acY'\^Y/'cT. 

Id 



A fiducial calculation is obtained from the parameters used previously for super- 
GZK estimates together with ac = (and E = 100 EeV). These values give 

Rmin = 0.26 Mpc. This minimum R corresponds to a maximum fossil density rjT = 
1/R^ = 60/ Mpc^ accumulated over T ^ 3 x 10^ yrs {cT = 100 Mpc). For the average 
galaxy density 0.01/Mpc^, this fossil density requires a burst frequency not more than 
one per 50,000 years per galaxy. (The minimum mean time between bursts per galaxy is 
< At O.OIR^T.) 

Note, however, the strong dependence of this minimum burst interval on the size N 
of the data set, it being proportional to R^ ~ 1/N'^. A five-fold increase in exposure 
would reduce the limit of 50,000 years down to 80 years. Detection of a discrete source 
would then be likely provided the mean time interval between bursts is not greater than 
80 years in a typical galaxy. Increasing exposure is especially advantageous in the search 
for discrete sources if the sources are isotropic bursts. 

Suh-GZK analysis: 

The background noise is not changed by the assumption of burst sources. It is still 
given by 

TV = VIujS = ViKrS. 

The signal-to-noise is therefore 

S/Af= , /y^=.[Il ^ 



ac 7]T \ K ac rjT r'^/^ 

For a specified value S for S/N, the maximum detection distance is r^. Using riT = 1/R^, 
r-s is given by 

rT = (24) 
The expected number of sources closer than this distance is 

= /d 47rr^ r]T dr = fu / AnaT]r^dr = T^foa'qr^, 
Jo Jo 
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which, upon substituting for becomes 



Detection is hkely, i.e. Prob{ro < r^) > 1/2, if /i > ln(2), which is 



R > i^y^' {cTf" {o^cf' {^f- (26) 



Using the estimate I£ = O.llTA^, this becomes v 



Rmm = 0.51(^)^/^(ac)V3(^)4/3. 

JD 



(27) 



A fiducial estimation can be done using ac = 1/{36E^) (for E = 10 EeV) and the 
same parameters that were used for steady isotropic sources: cT = 1500 Mpc, k = 
0.014 sr/Mpc, S = 3, = 10, 000, fo = 1/2. This gives Rmin = 1-3 Mpc. This minimum 
R corresponds to a maximum fossil density rjT = 1/R^ = O.S/Mpc^ accumulated over 
T ^ 4.5 X 10^ yrs {cT = 1500 Mpc). For the average galaxy density 0.01/Mpc^, this fossil 
density requires a burst frequency not more than one per 9.0 x 10^ years per galaxy. 

As in the super-GZK case, this minimum mean time between bursts in galaxies is 
proportional to 1/A^^. Modest increase in exposure can dramatically reduce the minimum 
time between bursts that would make the detection of a burst source likely. 



8 Jet bursts 

The previous two sections considered variations in which the cosmic ray sources emit in 
collimated jets or emit in isotropic bursts. It could also be that cosmic rays are emitted 
in collimated bursts. Both types of modifications to the basic model should then be 
incorporated together. The analysis here follows the now-familiar progression. 

Each burst emits a total number W of cosmic rays above some energy threshold, and 
they are emitted into a solid angle Qj. The spacetime density of bursts is t] bursts per 
unit volume per unit time. For an accumulation time T, the spatial density of bursts that 
contributed to the present cosmic ray intensity is n = rjT. Thus, 

Air Air I 

—I = nW = r]TW ^ W = . 

c c r]T 

The duration of a burst at distance r from the source is t = ar, and a detector with 
exposure S to that part of the sky will collect signal S if it is within the solid angle Q of 
the (magnetically spreading) jet, where 

^ _ W£ _ WE _ AttIS 
Qr'^t Qr'^ar ac rjT Qr^ 
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Super- GZK analysis: 



As in the earlier jet source analysis, suppose Q = Qj, i.e., the super-GZK particles 
retain their directions enough that the emission angle of the jet is approximately the solid 
angle of the jet's flux at the distance of Earth. Suppose K particles from the source are 
required for a signal-limited (super-GZK) detection. The distance is the distance at 
which a source is expected to produce that many signal events in a detector with exposure 
£: 

S = K ^ 4 = ^J^1^. (28) 

The expected number of sources within the volume of radius tk and exposed to the 
detector is 

^ = Id ! ^jr^V '^{r) dr = fo [ VLjaiqr^ dr = ^njar]r^ = ^VLjarf ( ^^J^ f'^. 

Jo io 4 4 ac rjT VtjK 

Using the mean separation of fossils that have contributed in time T (so = l/{riT)), 
the expected number /i is 

/^ = 4^(4vr)^/M7T^)^/Mfr/^4- ^29) 
4 \lj ac K CI 

As previously, the probability of the nearest source being detectable is greater than 1/2 
provided fj, > ln{2). For jet bursts, 

/. > lni2) ^ R> (i-)^/^ i^) (ac fi.)V3 (^)V3 (30) 
Making the approximation IS = 0.1 17A^ reduces this to 



Rrmn = ^-^{aC^jf^^Y'^cT. 

Id 



(31) 



A fiducial calculation can be done using the same adhoc values that were adopted previ- 
ously: cT = 100 Mpc, ac = 1/(36 x 100^), = 0.01 sr, K = 3, N = 50, fn = 1/2, 
giving Rmin = 0.024 Mpc. This is the mean separation of fossils of bursts that contributed 
to the cosmic ray density during the accumulation time T. That is, R = rjT. Using again 
the density of 0.01 galaxy/Mpc, the maximum mean time between bursts in each galaxy 
is given by < At >min= O.OIR^T, and for R = 0.024 Mpc, this gives < At >min= 41 yrs. 
For the fiducial parameter values adopted here, the nearest source is likely to be detectable 
provided the mean interval between bursts in each galaxy is at least 41 years. 

Sub- GZK analysis: 

As in the persistent jet sources, the solid angle of the jet is not approximately constant 
except at the highest energies. For the background-limited analysis, it increases with 
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distance (cf. Appendix C): Q = Sur. Substituting this into the expression for the jet 
burst signal gives 

^ ^ 47r IS 

3 ac rjT nr^ 

As usual, the background is luS, and the noise is its square root: TV = y/ InrE (using 
uo = Kr). The signal-to-noise ratio is therefore 



3 \ ac rjT r^/^ 

Setting a detection threshold S / M = gives the maximum distance at which this 
signal-to-noise ratio is expected: 

9/2 A-K ll£ 1 



S \ K,'* ac T]l h 

The mean density of jet bursts at distance r which are "on" and pointed at us is 
fo^TVL/ATi. Using T = ar and Vt = Snr, the expected number within a volume of radius 
r-s and exposed to the detector is 

= Id I o K rj a r ar = k a r] ry.. 

Jo 5 

Using the foregoing expression for r-£ and rjT = this becomes 

^ = ^(l![)10/9 (1)2/3 (1)1/9 1 (^^)5/9 (1)10/9 ^l/3_ (33) 

5 3 K ac cl L 

As usual, the nearest "on" source pointing at us is likely to be closer than and 
therefore detectable provided fi > ln{2). This condition is here equivalent to 

R > ih''^' i-^f" io^cf" {cTf T}''\ (34) 

3/d 47r I£ 

Substituting 0.117A^ for I£, this becomes 



Id 



(35) 



A fiducial calculation can be made again with the same parameter values as before: 
K = 0.014 sr/Mpc, N = 10,000, ac = 1/(36 x 10^), cT = 1500 Mpc, S = 3, /d = 1/2. 
The minimum jet burst fossil separation is needed for likely source detection is Rmin = 
1.3 Mpc. Using T ^ 4.5 x 10^ yrs, the minimum mean time between bursts in a single 
galaxy is < At >min= O.OIR^T = 1.0 x 10^ yrs. In the sub-GZK regime, detecting the 
nearest source is unlikely unless the mean time between bursts in a single galaxy is at 
least 100 million years. Note, however, that this calculated minimum time is extremely 
sensitive to the adopted value for cT, since it is proportional to R^T and R itself is 
proportional to (cT)^. A small change in the adopted parameter value for cT produces a 
dramatic change in the minimum mean time interval. 
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9 Luminosity function 



The basic model used in this paper assumes that cosmic ray sources are all the same in 
the sense that there is a single luminosity Q (particles emitted per unit time above the 
cosmic ray energy of interest). This is surely an overly simple idealization. There must 
be some distribution of luminosities n(> Q), i.e. the spatial density of cosmic ray sources 
with luminosity greater than Q. 

The luminosity function is constrained by the observed intensity /(> E) of cosmic 
rays. The luminosity function must account for the particle density — /: 



After imposing this normalization condition, a hypothetical luminosity function n(> Q) 
will generally not have a single parameter that could be used in place of R. The basic 
model in this paper has focused on the mean separation R between identical sources. 
There is no natural generalization of the R-analysis for an arbitrary luminosity function. 

Given any hypothetical luminosity function n{> Q), however, one can evaluate the 
probability P that one or more cosmic ray sources is detectable. It is given by P = 1 — e"'^, 
where /x is the expected number of detectable sources. Suppose a detector with exposure 
S requires K super-GZK events from some candidate source (or signal-to-noise E for sub- 
GZK events) to confirm a detection. At any distance r, there is a minimum Q{r) that is 
needed for that: 



for super-GZK or sub-GZK analysis, respectively. The expected number of detectable 
sources is then 



Jo 

which then gives the probability P = 1 — e~'^ that one or more sources is detectable for 
the hypothetical luminosity function n{> Q). 

10 Incomplete sky coverage 

If there are discrete sources of cosmic rays at all, they presumably surround us. The 
brightest nearby source could be in any part of the sky. It is important to achieve good 
exposure to the full celestial sphere. No cosmic ray observatory has so far been built with 
full-sky exposure, although the Auger Observatory has been designed for that. 

The analyses in this paper have used a simplifying approximation that a single-site 
detector has the same good exposure to a fraction of the sky denoted by fo- The southern 
site of the Auger Observatory, by itself, has good exposure to approximately one quarter 
of the sky at the most southern declinations, decreasing exposure over half of the sky, 
and zero exposure to the northernmost quarter of the sky. Adding the complementary 
exposure from its northern site would yield nearly uniform acceptance to cosmic rays from 
all parts of the sky [9] . 





> K or 
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The formulas in this paper become more comphcated for the real situation in which 
exposure S{u) is not constant over the celestial sphere but depends on direction u. (Here 
u denotes a unit direction vector.) In the case of isotropic persistent sources, for example, 
the signal is 

IR^ 

(See equation [21) For the super-GZK analysis, the expected signal is K at radius given 
by 

, IR^ 



^ cTK 



£(u). 



The expected number of detectable sources out to distance tk, which now depends on 
direction u, is 



/i 








"^"^^r^rfrrffi=l(^)^/^ l{8{u)f'Hn. (36) 



The nearest source is still more likely to be detectable than not provided /i > ln{2), which 
now becomes 

3ln2 2/3 (^TK 
^ ^ J{S{u)Y/^dQ ' I ■ 

Similar modifications pertain to the sub-GZK analysis. Both the signal and the noise 
depend on direction u, and the ratio is 



The distance at which this has a prescribed value S is given by 



r 



/ R^ 



5/2 



^fW)■ 



The expected number of detectable sources is 

^ = / /"^"^ y'<ir = y-^r/' {^f' J {S{u)f''dn. (37) 
The condition /i > ln{2) becomes 

R > ,3 lu2f'^ ^ (cT^? (7^£^)^'= 

A simple schematic approximation for the Auger South exposure is £^ = Sog{v) where 
V = sin{dedination) , Sq is the rich exposure in the southern sky, and g{v) is the simple 
function 

1 if -1 < < -1/2 

giv)={ 1/2 -V if -1/2 < i; < 1/2 
ift;>l/2. 
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(Note that g{v) + g{—v) = 1, so identical sites in the north and south provide uniform 
sky coverage Sq in this approximation.) 

This schematic model of exposure for Auger South allows an approximate analytic 
evaluation of J{S{u)ydQ for the various powers 7 which would arise in the different 
models considered in this paper. The resulting formulas are certainly less transparent, 
however, than for the ideal detector with uniform exposure. 

An effective sky coverage fr, for Auger South can be defined by requiring that the 
expression for /j, in equation [36] be a scaled version of the expression in equation HI so 



For the sub-GZK case, the analogous condition on coming from equations [37| and [8] 
leads to 



These results suggest that 1/2 is a suitable estimate for the effective sky coverage fn that 
pertains to a single-site observatory. 

11 Summary and conclusions 

The likelihood of detecting a discrete source of ultra-high energy cosmic rays depends 
on many variables with unknown values. It has here been assumed that cosmic rays are 
produced in sources that are randomly distributed throughout the universe with some 
mean separation R. The analyses have focused on the question, "What condition on 
the mean separation ensures that the detection of the nearest source has more than a 
50% chance of being detectable?" The question is simplified by assuming that all sources 
have the same luminosity Q. The answer certainly depends on the amount of exposure S 
that the detector has to the sources. Since the number of arrival directions in a data set 
increases in proportion to the detector's celestial exposure, the answer can be regarded 
alternatively as depending on the number N of arrival directions collected above the 
chosen energy cut. It also depends on the cosmic ray accumulation time T for that energy 
cut. Another variable that affects the answer is the statistical significance required for 
detection - the required number of events K in the signal-limited case (super-GZK) or 
the required signal-to-noise ratio S = S/Af in the background-limited case (sub-GZK). 
In addition, the answer depends on whether the sources are permanent sources emitting 
isotropically, permanent sources with beamed emission, isotropic bursts, or jet bursts. 
Magnetic fields also have an important impact on the answer for the background-limited 
cases. 

Properties of intergalactic magnetic fields are not well established. For permanent 
sources, the relevant information about magnetic fields is summarized by the coefficient 
K which governs how the (Gaussian) solid angle of arrival directions grows with distance 



An IEqR 
Y^cTK 




This yields 
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from a source: uj = nr. Appendix B examines how k can be calculated from the mean 
field strength and coherence length in a simple model. 

For transient (burst) sources, magnetic smearing of arrival times is relevant as well as 
the smearing of arrival directions. The time smearing is encapsulated in the coefficient a 
hy T = ar. Appendix C shows how to estimate a from the magnetic field properties in 
the simple model that is used to estimate k. 

The tables below collect the boxed formulas appearing in the text. For each class 
of models, there are two formulas for the minimum mean separation Rmin such that the 
probability of detection is greater than 1/2. One formula pertains to the super-GZK 
(signal-limited) regime, and the other formula pertains to the sub-GZK (background- 
limited) regime. 

Table 1: Formulas for Rmin- This is the minimum mean source separation for which 
one expects the nearest source to be detectable. Formulas are tabulated for the signal- 
limited (super-GZK) case and the background-limited (sub-GZK) cases. These formulas 
are highlighted by boxes in the text. For a full-sky observatory with uniform exposure, 
should be half the number of arrival directions in the data set of events above the energy 
corresponding to the cosmic ray accumulation time T, and the fraction of sky exposed to 
the detector is //) = 1. For a single site like Auger South, A^ is the total number of events 
in the data set, and fo ~ 1/2. For non-uniform exposure, it would be better to replace A^ 
by /£'/0.117 and use the exposure £ that pertains to the target source celestial position 
and the chosen energy cut for which the cosmic ray intensity is /. 



Source Type 


Super-GZK 


Sub-GZK 


Steady isotropic 


R ■ = 4 1 — cT 


Rm^n = 1.4 (cT S)2 f 


Steady Jets 




Rmm = 1.6(^)^/3(cT)8/3«:5/3 


Isotropic Bursts 


Rmin fi3 ('^'-') ^ (iv^ ^ 


/2_ = 0.51(f)^/3(ac)V3(^)4/3 


Jet Bursts 


Rm^n = g (acfi j) ^3 (|)4/3cT 


Rm^n = '-f («c)V3 (f )5/3 ^2 (^y)3 



A statistical detection of multiple sources can be expected prior to the detection of 
any one discrete source. Study of the two-point correlation function (or, equivalently the 
angular power spectrum) for ultra-high-energy cosmic ray arrival directions can exhibit 
evidence for many poor clusters of arrival directions even if there is not any one rich cluster 
of arrival directions that is individually detectable [TU]. Evidence for an autocorrelation 
in the AGASA data was published [11]. Moreover, a correlation of arrival directions with 
a catalog of candidate sources can supply evidence of discrete sources even if there is 
no statistical evidence for clustering of the cosmic rays themselves. Exploratory searches 
have also produced evidence for correlations of that type [12]. These specific claims will 
be thoroughly tested using larger data sets obtained with new observatories. The failure 
to detect any individual discrete source so far suggests that there is no really bright cosmic 
ray source in the sky. It is therefore likely that sources will show up collectively in one of 
these ways before any individual source becomes obvious. 

Whether or not the brightest cosmic ray source will soon become detectable depends 
on many unknown parameter values. The issue must be decided observationally. Fiducial 
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estimates show that it is a close calL The answer can go either way, depending on 
assumptions about the unknown parameter values and the nature of the sources (e.g. 
bursts, jets). Present detector exposures are already obtaining important constraints on 
the unknown astrophysical parameters. The viable parameter space will shrink rapidly 
as the exposure increases. The search will be especially rewarding, however, if it leads to 
the study of one or more sources and the intervening magnetic fields. 
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13 Appendix A: Minimum and maximum R 

The analysis in the text focuses on the mean separation R between identical sources 
(density of sources n = 1/R?). The nearest source is unlikely to be detectable if R is too 
small, because then each source is too weak. For super-GZK analysis, there is an upper 
limit Rmax as well as lower limit Rmm stemming from the fact that the nearest source 
is unlikely to be much closer to us than the mean separation i?, so R cannot be much 
greater than cT, where cT is the effective survival distance for those super-GZK particles. 
In a volume of radius r around us, the number of detectable sources is 

, = '-^r^/R\ (38) 

As explained in section 3, the nearest source is likely to be detectable provided /i > ln2. 
The critical value /i = ln2 gives Rmin and Rmax- Equation [38] shows that these are simply 
proportional to the corresponding volume radii: 

Rmin Xv^iji and Rmax •^'"max ; 



where the proportionality constant is 



47r/DM/3 



Equation [T] gives the luminosity Q = ^^^^^ per source with mean separation R which 
collectively account for the cosmic ray density Airl/c. The signal S seen at distance r 
after exposure £ in the basic (isotropic) model is 

cTr'^ ' 

where the GZK attenuation factor e~^^^'^ is here included. If K events are required for 
a source detection, the critical case is obtained by setting S = K. Measure distances 
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relative to cT, so r = ^cT and R = X^cT when the distance r is = Tmm or = Vmax- 
Then the equation above for S reduces to 

f.-^ = ^. (39) 

A solution for requires -^^^ < i, and there are two solutions except if the equality holds. 
One solution has ^ < 1 (i.e. vk < cT), and the other has ^ > 1 (i.e. rx > cT). For cases 
in which <^ 1/e, then = e"''^/'^-^ = 1. For r < r^, the attenuation factor e~'^^^'^'^ 
is then very close to 1, and the simpler analysis following equation [2] is fully justified. 

The approximation IS = 0.1 17A^ introduced in section 3 can be used here. Adopting 
/z) = 1/2 as in the fiducial calculations, one gets 

K 

<l/e^K < 0.13N. 



Therefore, provided the data set has more than 8 times the number of arrival directions 
K needed for detection of the discrete source, the analysis in the text (omitting the e"^/"^"^ 
attenuation factor) is adequate. 

The upper limit Rmax = X^cT is given by the other solution of equation [39] for which 



14 Appendix B: Blurring by magnetic fields 

The Larmor radius p characterizes the bending of charged particle trajectories by magnetic 
fields. For relativistic particles, it is given by 



eV 



PMpc 



= for extragalactic applications; 



300 Z Bq 1 p^^^ j_ ^Eev fQj, galactic applications. 

The approximate equations on the right follow from the exact equation on the left using 
the relations: 1 EeV = 10^^ eV, 1 nG = IQ-^ G, 1 Mpc=3.1 x 10^^ cm, 1 = 
10~^ G, 1 kpc=3.1 X 10^^ cm. The particle's electric charge Z is in units of 1 proton 
charge. In traveling distance D through a perpendicular B-field, a particle's trajectory is 
bent by the angle 6 = D/p in radians. 

The Galaxy has a regular magnetic field which tends to follow the spiral arms, and 
also superposed irregular fields which change direction over short distances along any 
path. Particles arriving from a distant point source will be systematically deflected by 
the regular field. An estimate is that they will encounter, on average, a perpendicular 
magnetic field of about 2 pG acting over a path of roughly 1 kpc. Trajectories from a 
single source would then be systematically bent through the angle 

, I kpc (lkpc)Z(2nG) 
6 = = radians. 



P E, 



EeV 
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For protons {Z = 1) this is about 12° at 10 EeV and 1.2° at 100 EeV. 

A uniform field acting along the entire path from source to detector would cause the 
arrival direction to differ from the source direction by only 6/2. This is because, in that 
case, the particle does not start from the source in our direction; its initial direction also 
differs from our line of sight by half of the trajectory bending angle {0/2). For a distant 
extragalactic source, however, the regular field changes the arrival direction by the full 
angle 9 from the particle's direction of entry into the galaxy (which is the direction from 
detector to source if the source is very distant and the particle travels on a straight line 
until reaching the Galaxy). 

All protons of one energy are deflected the same amount by the regular magnetic field. 
With arrival directions of two or more protons of measured energies from the same source, 
one can determine the product BD and the source direction. Here BD is the transverse 
magnetic field integrated over its range along the incoming trajectory. The regular field 
produces an arc of arrival directions on the sky ending at the source direction [E = oo) 
with E'-dependent angular offsets of ^ = BD/ E. 

Because of the steep energy spectrum, it can be expected that roughly 3/4 of the 
arriving particles will have energy less than twice the analysis threshold energy. Using the 
above estimate of deflection by the regular galactic magnetic field, one would expect that 
above 10 EeV, for example, the regular field should spread 75% of the arrival directions 
along an arc of roughly 6 degrees (plus or minus 3 degrees), with the center of that arc 
displaced from the source direction by approximately 9 degrees. Above 100 EeV, the arc 
would be roughly 0.6 degree (plus or minus 0.3 degree) with its center displaced from 
the source by approximately 0.9 degree. The input values for these estimates (2 /iG for 
the transverse regular field strength and 1 kpc for the effective path length) are crude 
estimates, and actual values depend critically on the direction to any given source. This 
rough estimated does indicate, however, that clusters of arrival directions should not be 
destroyed by the Galaxy's regular magnetic field for most source directions. 

Irregular magnetic fields are the other concern in charged particle astronomy. These 
are fields that do not have a consistent direction over any particle's trajectory. The 
particle's direction is continuously being deflected by small magnetic bends that cause it 
to meander in a random- walk manner. The result is formally the same as multiple coulomb 
scattering of energetic charged particles in matter. A uni-directional initial beam becomes 
a Gaussian distribution of particle directions centered on the undeflected direction. The 
width of the Gaussian distribution (a) increases in proportion to the square root of the 
path length r. Appendix C shows that the effective solid angle is (a; = Aira"^), which 
increases in proportion to the path length. The proportionality constant k, is defined by 
u = nr. 

A simple model of magnetic blurring by irregular fields is that a particle encounters 
a different field orientation in each segment of length L along its path. The field has a 
mean strength B with a random orientation which is constant over each segment. Consider 
deflection in any plane containing the original direction. The B-component perpendicular 
to that plane is expected to be B/^%. Over a segment length L the deflection in the 
plane is 
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This is the random walk step size. After n := r / L steps, the directions in that plane are 
distributed with a Gaussian of width 



[7 LZB rL ZB 

This is the distribution of directions relative to an original beam direction. What is 
relevant is the direction of a detected particle relative to the direction from the detector 
to the source. That direction to the source does not correspond to the original particle 
direction because the particle will have been displaced laterally from the beam. The lateral 
displacement is correlated with the offset of the arrival direction. As in multiple Coulomb 
scattering, the arrival directions are distributed around the direction to the beam origin 
with a Gaussian a which is smaller by 1/a/3. This is a straightforward consequence of 
the statistical correlation of angular offset with spatial offset in the random walk process. 
The direction back to the source is different for particles that have a net deflection to the 
left than for for those that have a net deflection to the right, and the final direction is 
statistically correlated with the net deflection. 

Using K = uj/r together with uj = Ana'^ and a = ^/j-^ yields a formula for k: 

An LZ^B^ 

K = . 

9 

This expression for k is used in the text for evaluating the magnetic blurring due to 
random intergalactic magnetic fields between a distant source and the Earth. There is 
also a contribution by irregular magnetic fields within the Galaxy, and this expression for 
K can be used to estimate its magnitude. Suppose B = 3 fiG for the randomly-oriented 
field strength, that L = 0.1 kpc, and the travel path through those irregular galactic fields 
has length r = 1 kpc. Then the Gaussian spread a is given (for Z = 1) in radians by 



a = a/ Ld/An = \j 
Setting E = 10 EeV gives a = 1.8°. 



LZ^B^r _ 0.32 
9^2 ~ 



15 Appendix C: Time spreading by magnetic fields 

The simple model of irregular fields used in Appendix B provides an estimate for the 
coefficient a that governs the time spread r for a source at distance r by r = ar. The 
estimate for r is based on the expected difference in transit time for a charged particle 
of energy E compared to an undeflected neutral speed-of-light particle. There is a time 
difference in each path segment of length L = p6 because the curved trajectory is longer 
than the straight line distance between the endpoints. (Here p is the Larmor radius based 
on the perpendicular magnetic field, and theta is the trajectory bending angle while 
traveling distance L.) The time difference is 

At = -{L - 2psin{9/2)) = ^{9 - 2sin{9/2)) ^ ^— 



c c c24 24cp2 
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The total time difference is then 




Using R = E/ZBx and < Bl 



>= |-B^, this gives the coefficient a = r/r: 



a = 



36cE2 ■ 



Here B is the magnetic field strength whose direction is randomly oriented along each 
trajectory segment of length L. 

This estimate ignores a second-order contribution to r due to the segments themselves 
meandering about the straight line from the source to the arrival point. 

16 Appendix D: Signal and noise with Gaussian-distributed 
arrival directions 

Random intergalactic magnetic fields are expected to produce a Gaussian distribution of 
arrival directions at the detector, centered on the source direction. Suppose an observed 
source has arrival directions distributed about a central direction with a Gaussian of 
width a. This means in any one dimension the probability distribution for offset 6^ is 
P{9x) = y=^exp(— 6*^/2(7^), and the 2-dimensional (space angle) offset 9 = ^/O^ + 9y has 
probability P{0) = ^;^exp{—9'^/2a^). When testing a discrete source with an apparent 
Gaussian distribution of arrival directions, it is sensible to give more weight to arrival 
directions that are near the center of the distribution and little weight to arrival directions 
that are far from it. The appropriate weighting function is 



This choice enjoys three important properties: 

(1) The shape of the Gaussian function of width a maximizes the signal among all possible 
weighting functions of the same integral normalization. 

(2) With the Ana"^ integral normalization, the weighted integral of any uniform back- 
ground has expected fluctuations given by the square root of the background itself. This 
expected background is Ana^ x {background density), and the RMS fluctuation in that 
weighted background integral is its square root. 

(3) With the Ana"^ integral normalization, a density NqP{6) of smeared-out arrival direc- 
tions due to A^'o cosmic rays from a discrete source results in a weighted integral equal to 
A'^o, i.e. the actual number of smeared-out directions. 

Defining signal and background as weighted integrals with this weighting function, 
their difference is the expected number of events producing the signal. Moreover, the 



w = A7ra^P{9) = 2exp{-9^ /2a^). 
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usual S/N statistical significance pertains with the noise A/" being simply the square root 
of the background, as in analyses without a weighting function. One can regard 2a as an 
effective radius, giving Aira"^ as an effective collecting area. 

A more careful analysis should use Fisher distributions [l3] rather than Gaussians for 
celestial analyses, especially for the broad distributions that are expected in sub-GZK 
analyses. 
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